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SUMMARY 


A double-acting electromagnetic thrust bearing is normally used to counter the axial loads in many 
rotating machines that employ magnetic bearings. It essentially consists of an actuator and drive 
electronics. Existing thrust bearing design programs are based on several assumptions. These 
assumptions, however, are often violated in practice. For example, no distinction is made between 
maximum external loads and maximum bearing forces, which are assumed to be identical. Furthermore, 
it is assumed that the maximum flux density in the air gap occurs at the nominal gap position of the 
thrust runner. The purpose of this paper is to present a clear theoretical basis for the design of the 
electromagnetic thrust bearing which obviates such assumptions. 


INTRODUCTION 


The basic design analysis of an electromagnetic thrust bearing is well known (1 , 2, 3, 4). In these 
analyses, the maximum bearing force is generally assumed to equal the external load. However, this 
assumption ignores the inertia forces due to the vibrating shaft. These inertial forces may be larger than 
the external loads under certain situations, thus invalidating the conventional design analysis. 

Another assumption normally used is that the maximum force requirement occurs at the nominal or 
equal gap position, when the thrust runner is centered between the two stators. This assumption holds 
only if the maximum movement of the runner is negligibly small compared to the air gap. An adequate 
design should take the motion of the runner into account when computing the maximum ampere turns 
required. 


NOMENCLATURE 


A p 

single pole face area [m 2 ] 

I slut 

static component of coil current [A] 

B u B 2 

flux density in left and right 

J 

maximum coil current density 


electromagnet respectively [T] 


[A/m 2 ] 

C br$ 

bearing damping coefficient [N-s/m] 

Is 

axial length of slot [m] 

kbrg 

bearing stiffness coefficient [N/m] 

N 

number of turns in coil 

m 

total moving mass (shaft and runner) 
[kg] 

R 

resistance of coil [fi] 



t 

time [s] 

f 

coil copper factor 

to 

instant when bearing force is 

fp 

coil packing factor 


maximum [s] 

I'U Fi 

force on runner from left and right 
electromagnet respectively [N] 

V 

power supply voltage [V] 



X 

axial displacement of runner [m] 


net bearing force on runner [N] 





amplitude of axial displacement of 

ciyn 

amplitude of dynamic component of 
external force on shaft [N] 


runner [m] 



a 

correction factor accounting for flux 

/' ext 

external force on shaft [N] 


leakage and fringing 

^ s tat 

static component of external force on 
shaft [N] 

mean air gap [m] 

4 

phase lag of runner displacement 
with respect to dynamic component 
of external force [rad] 

#0 

o 

air gap for left and right 

Flux in magnetic circuit [Weber] 

gi.ft 


hs 

electromagnet respectively [m] 
radial height of slot [m] 

fto 

permeability of free space [4n x 10‘ 7 
N/A 2 ] 

h,h 

current in the left and right 
electromagnet coil respectively [A] 

0) 

frequency of dynamic external 
force [rad/s] 

h 

bias current in coil [A] 


phase lag of coil current with respect 
to dynamic component of external 

Idyn 

amplitude of the dynamic component 
of coil current [A] 


force [rad] 


I max maximum coil current [A] 
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THEORY 


A double-acting electromagnetic thrust bearing, Fig. 1, is made up of two electromagnets (stators), 
one on each side of the thrust runner (rotor) and separated by an air gap. Such a system is inherently 
unstable. It is stabilized by sensing axial position and using this information to control the current in the 
stators via an electronic controller. The controller is typically of the proportional-integral-derivative 
(PID) type, and may be implemented in analog or digital form. 


Force Capability 


Each of the stators applies an attractive force on the runner. The free-body diagram of the thrust 
runner is shown in Fig. 2. For dynamic equilibrium, the equation of motion is 


mx = F ext - F brg 


( 1 ) 


where the external force and the resulting motion are 

( 2 ) 
(3) 


F eXt = Kat + Fjyn Soffit 
X = x dyn sin (<ot - <£) 


On substitution, we get 

F hrg = F ext - mx = F stai + F dyn sin ojt + mco 2 x dyn sin(ra t - <£) (4) 


The bearing has to be designed to accommodate the vector sum of the external force and the inertia 
force. The phase lag determines whether the inertia force adds to or subtracts from the external dynamic 
force. Fig. 3. For subcritical operation {</> < k/2), neglecting the inertia force can result in an undersized 
thrust bearing that “bottoms out” during operation. For supercritical operation ( <f> > n/2 ), the inertia 
force fights the external dynamic force, and ignoring this can result in an overdesigned bearing. If the 
phase lag is small 


(f> « k/2 ( 5 ) 

the bearing must be designed for a load capacity of 

= F sta, + F dyn + ^ X dyn ( 6 ) 

This force requirement occurs at cot = n/2 . Since the two electromagnets act in pull-pull mode, F, and 
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f' 2 are always positive. The net bearing force is the difference of the two: 


F brg = F, - F 2 (7) 

The bearing force is at a maximum when one electromagnet is applying its maximum pull force and the 
other is turned off. Thus, for the double-acting arrangement, we must have 


l,max 


= F. 


2, max 


= F 


brg, max 


( 8 ) 


Each electromagnet must be designed for this load capability. For the special case where F exl =0, the 
bearing must still be capable of withstanding the inertia force at the natural frequency of axial vibration 

of the system. System design is then dictated by desired transient response, characterized by overshoot 
and settling time (5). 


Ampere Turns 


The thrust bearing consists of two electromagnetically biased and excited magnetic circuits with two 
air gaps per circuit. Each circuit has an outer and an inner pole in the stator, and back iron to complete 
the flux path in the stator. Fig. 1 The flux lines traverse the air gaps and complete their path in the 
runner. Both the runner and the stator are made of magnetically permeable material. The pole face areas 
are typically made equal in order to ensure uniform flux density in the stator. 

The bearing force is related to the flux densities in the two electromagnets by 

(9) 

while the two air gaps are given by 


Sx = go +X 

g% ~ go — X 

The corresponding currents in the two electromagnets are 

I \=K+Kuu + 1 dm sin(a# - (//) 
K = K ~ Km, ~ Sin (^ ~ ?) 


( 10 ) 


( 11 ) 


The maximum external force, F mmax occurs when cot = tc'2. This is given by 

F — F I _ 17 . i/ 

ext, max ext\ cot = 7Z , 2 1 stat ^ r dyn 


(12) 
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However, F/, rg may not reach its maximum at the same instant due to the effect of rotor inertia. Let the 
bearing force reach a maximum at some time t 0 . At this time / 0 , then, we must also have, from (9) 





(13) 


Since 


B 2 > 0, B\> 0 (14) 

the difference reaches a maximum only when B 2 - 0 and B\ is at a maximum, or B\ = 0 and Bi is at a 
maximum. Let’s consider the former case: 


B 2 = 0 => 1 2 = 0 => I stat +1^, sin(a>T -y/)= / h = / max /2 

/,=/* + l s, a , + U. sm(«7'- y/) = I h + I h = / max 


(VoAO 

|/« 


( Mo N^ / max 

V 2a J 

Mi 

t-T 

l 2a J g 0 + x dyn sin (coT - <f) 


(15) 

(16) 

(17) 

(18) 


We can rearrange ( 1 8) to determine the (NI) m ax required to produce the magnetic flux density, B 1 max : 

( ■ NI ) max = a [g 0 + jy, sin(<y7’ - <*)] (19) 

Bo 


When the material operates up to saturation, B max B sal , and (NI) max can be calculated, provided tO is 
known. The correction factor a may be determined using the techniques discussed in (2). The pole face 
area required may be determined using Eqn. (15) above. 

The time / 0 when the bearing force developed is maximum can be determined analytically, using 
Eqn. (4). It is given by 


T - (x-S)/to 


(20) 


Here, the angle S can be calculated from 


S = 


tan 


" Ty, +rna) 1 x^ l cos£ 

v mc ° 2 x dyn sm</> J 


( 21 ) 


When the small phase lag condition of Eqn. (5) holds, the maximum force occurs at £ * x dyn . This 
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implies that the designer must ensure sufficient ampere turns to saturate the magnetic material of the left 
electromagnet for its maximum gap position of the runner. In general, however, the phase angle <t> is not 
known beforehand, since the bearing design affects the system dynamics. Thus, an iterative process 

must be employed or bearing stiffness and damping parameters assumed in order to determine this 
phase: 

F br g = c brg x + k brg x (22) 

For a given dynamic to static load ratio and a required vibration to gap ratio, the required stiffness may 
be calculated using the methodology outlined in [6], 


Number of Turns and Maximum Current 


Assuming a supply voltage v driving the coil of the left electromagnet, we must have 


V= N^- + I [ R 
dt 


( 23 ) 


The resistive load is typically small compared to the inductive load, and may be neglected. Then, 


<t> 


= — f vdt = — V si ncot 

N J coN 


( 24 ) 


assuming a sinusoidal supply voltage of the form v = V mn cos cot [7], The maximum flux level is 
therefore 


O = B A - 

max max p 


V. 


coN 


( 25 ) 


The power supply voltage required is often determined by force slew rate requirements [8]. However, 
y max is limited by available power supplies that operate at the frequency of interest. The number of turns 
is thus fixed as 


N = 


V 

max 
°> A p B m ax 


This, in conjunction with Eqn. (19), fixes the maximum coil current: 

(NI) 

V / max 


^ ma\ 


N 


( 26 ) 


( 27 ) 
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The maximum coil current density is typically fixed at around J = 5 x 10 6 A/m 2 , using which the coil 
wire gage may be obtained 


A., 


u 


(28) 


The coil copper factor f. — the ratio of bare copper cross-section to wire cross-section — is typically 
less than 0.7. The coil occupies a slot of axial length /, and radial height h„ so that the slot cross-section 
is 


A=l.h . 


(NI) 

V / mai 

1 7 ? 


(29) 


EXAMPLE — HIGH SPEED FLYWHEEL 


Let us consider a flywheel of mass 1 kg (2.2 lb) rotating at 6283 rad/s (60,000 RPM) and supported 
on magnetic thrust bearings. Let the magnitude of the dynamic external force be 445 N (100 lbf), so that 

F ext = 445 sin 6283i 

It is desired that the maximum excursion of the flywheel under this loading be limited to 2.5 x 10 ' m 
(1 mils). So, the runner motion is 


x = 5 x 10 6 sin(6283£ - 


If <|> satisfies Eqn. (5), we must design for 

If, max I = If max I = |F, max I = 445 N + 1974 N = 2419 N 

I l.max I max I I or^,max I 

as required by Eqn. (6). It is evident that a bearing designed to withstand only the 445 N of dynamic 
external force would be inadequate for this system. 

If the mean air gap is go = 2.5 x 10 -4 m (10 mils), B max = 1.0Tand Eqn. (19) is applied now, the 
maximum ampere turns required is 


(NI) - 438 A - turns 

assuming that the correction factor is a = 1 . This is about 40 A-tums more than for conventional 
designs, which assume saturation at the equal gap position. Also, using Eqn. (15), the pole-face area 
required is 
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A p = 5.59 x 1(T 4 m 2 (0.87in 2 ) 


We can now determine the number of turns N by applying Eqn. (26) and assuming a supply voltage 
of, say, 200 V: 


N = 57 turns 


The maximum current in the coil is 


I 


max 


7.7 A 


CONCLUSION 


We have shown that two assumptions often employed in designing electromagnetic thrust bearings 
can lead to inadequate products. The inertia force must be taken into account when sizing the thrust 
bearing for load capacity. For subcritical operation, this inertia force increases the bearing load capacity 
requirement. However, if one designs for supercritical operation, the inertia force fights the external 
force and reduces the load capacity requirement. The bearing size may be reduced by taking advantage 
of this fact. 

Moreover, during the vibration of the thrust runner in response to an excitation force, the maximum 
bearing load capacity is required at a position that is different from the nominal gap position. The 
implication is that the ampere turns required is more than that for the maximum-force-at-the-nominal- 
gap assumption. The maximum available voltage supply to energize the electromagnets then determines 
the maximum coil current and the number of turns. 
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Figure 1 Double-Acting Electromagnetic Thrust Bearing 
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Figure 2. Free-Body Diagram of Thrust Runner for a 
Magnetic Thrust Bearing 
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F dyn 



(a) Subcritical Operation 


^dyn 



(b) Supercritical Operation 


Figure 3. Phasor Diagrams for Motion of Thrust Runner 
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